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ABSTRACT 


We  have  made  a  semi -empirical  analysis  of  our  extensive 
body  of  observed  cosmic  ray  intensity  data  from  Pioneers  10  and 

cj  £K  <  r  f  r  r 

11  and  related  spectral  information  from  other  .author  in  order 
to  infer  the  radius  R  of  the  modulation  region  surrounding  the 
sun.  During  the  period  1972-1985  the  inferred  values  of  R  vary 
with  time  systematically  and  in  a  manner  generally  similar  to 

that  of  sunspot  numbers.  The  range  of  values  of  R  is  from  42  AU 

obe  ><  T 

at  the  time  of  minimum  solar  activity  (~  1976)  to  88  AU  about 
1.5  yr  following  the  time  of  maximum  solar  activity  (  -v  1980). 


A  specific,  testable  prediction  is  that  Pioneer  10  will  cross  the 
modulation  boundary  in  1988  and  will  remain  in  its  vicinity  for 
several  years  thereafter 
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One  of  the  classical  problems  of  heliospheric  physics  is 
the  observational  determination  of  the  boundary  of  the  region 
around  the  sun  within  which  the  11-year  periodic  modulation  of 
galactic  cosmic  ray  intensity  is  caused  by  the  entrained  magnetic 
field  in  the  outflowing  solar  wind.  This  boundary  has  not  yet 
been  encountered  by  Pioneer  10  at  its  present  heliocentric  dis¬ 
tance  r  of  38  AU  in  the  solar  antipex  direction  nor  by  Pioneer  11 
at  21  AU  in  the  solar  apex  direction. 

The  scale  size  of  the  region  is  the  reciprocal  of  the  ob¬ 
served  mean  radial  gradient  of  intensity,  namely  50  AU.  Crude 
physical  arguments  suggest  that  the  radius  R  of  the  modulation 
boundary  lies  between  50  and  100  AU  [Van  Allen  and  Randall, 

1985],  Uebber  and  Lockwood  [1981]  estimate  that  R  >  65  AU;  and 
McKibben  et  al.  [1982]  estimate  that  R  lies  in  the  range  50  to 
70  AU  at  both  solar  maximum  and  solar  minimum.  Both  groups  of 
authors  use  observed  interplanetary  data,  modulation  theory,  and 
an  assumed  interstellar  energy  spectrum  as  the  basis  for  their 
estimates.  Kurth  et  al.  [1984]  have  interpreted  a  distinctive 
episode  (1983-84)  of  interplanetary  radio  emission  in  the  fre¬ 
quency  range  2  to  3  kHz  as  suggesting  that  the  heliopause 


(not  accurately  identical  to  the  modulation  boundary)  lies  at 
R  »  46  AU. 

The  new  aspect  of  the  problem  that  is  exploited  herein  is 
the  recent  availability  of  an  extensive  body  of  cosmic  ray 
intensity  data  (Ep  >  80  MeV/nucleon)  from  Pioneers  10  and  11  over 
a  complete  cycle  of  solar  activity  and  to  great  distances  from 
the  sun  [Van  Allen  and  Randall,  1985].  The  intensity  is  a 
function  of  both  radial  distance  r  and  phase  of  the  solar 
activity  cycle  (as  well  as  of  sporadic  solar  activity)  (Figure  1) 
Gleeson  and  Axford  [1967,  1968]  have  shown  that  the  convection- 
diffusion  model  of  cosmic  ray  modulation  [Parker,  1963,  1965] 
yields,  under  certain  simplifying  assumptions,  a  one-parameter 
characterization  of  the  relationship  of  the  differential 
intensity  of  a  particular  species  of  particle  at  infinity  (i.e., 
in  the  interstellar  medium)  to  that  at  an  observing  point 
interior  to  the  modulation  boundary  [cf.  Ehmert,  1960,  and  Freier 
and  Uaddington,  1965] . 

We  use  this  approximate  formulation  together  with  our 
integral  intensity  data  and  spectral  data  from  other  authors  to 
make  speculative  estimates  of  the  heliocentric  radius  of  the 
boundary  of  the  modulation  region  as  a  function  of  phase  of  the 
solar  cycle. 
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DEMODULATION  MODEL 

According  to  the  formalism  of  the  Gleeson-Axford  theory 
[see  also  Goldstein  et  al.,  1970],  the  effective  change  in  energy 
of  a  charged  particle  that  moves  from  a  source  point  P  on  the 
boundary  along  a  dynamical  trajectory  to  an  interior  point  Q  is 
-Ze4,  irrespective  of  its  path  from  P  to  Q;  that  is,  a  particle 
whose  observed  kinetic  energy  at  Q  is  T'  must  have  had  a  kinetic 
energy  (T'  +  Ze4)  at  P.  The  unidirectional  differential  (in 
energy)  intensity  J'q(T')  at  Q  along  a  dynamical  trajectory 
leading  from  P  to  Q  is  related  to  that  at  P  (in  a  quasi -steady 
state)  by 


Oq(T’) 


ip(T'  +  ZeO 


T’(T’  +  2Amc  )  (T’  +  Ze  * )  (T’  +  Ze  *  +  2Amc  ) 


(1) 


where  A  is  the  atomic  mass  number  of  a  particular  species  of 
energetic  ion,  Ze  is  its  charge,  me  is  the  rest  energy  of  the 
proton,  T'  is  the  kinetic  energy  of  the  ion  at  the  observing 
point  Q  and  (T'  +  Ze  4)  is  its  kinetic  energy  at  the  boundary. 
(Note  that  4  is  positive  for  ions  and  negative  for  electrons.) 


Denote  the  energy  per  nucleon  by  T  such  that  T'  -  AT.  Then 


(1)  becomes 


VT>  .  VT  *  jr]_ 

T(T  +  2mc2)  =  (T  ♦  )  (T  ♦  ~  ♦  ) 


(2) 


In  the  present  study  we  use  observed  data  only  for  protons 
(Z/A  -  1)  and  alpha  particles  (Z/A  -  0.5). 


The  high  energy  (T  >  30,000  MeV/nucleon)  portions  of  the 
differential  spectra  of  primary  protons  and  alpha  particles  at 


1  AU  [Meyer,  1969]  [Ryan  et  al.,  1972]  [Fulks,  1975]  can  be  well 
represented  by: 


with  spectral  index  7  =»  2.7. 

We  assume  that  the  spectral  form  of  (3)  is  valid  down  to 
T  •»  80  MeV/nucleon  in  the  nearby  Interstellar  medium.  This 
assumption  is  plausible  and  there  is  no  present  evidence  against 
it  but  it  must  be  regarded  as  tentative  and  may  later  be  found  to 
be  wrong. 

Then  combining  (2)  and  (3)  and  denoting  Ze$/A  by  <p, 
we  have  the  following  two  working  formulae  for  protons  and  alpha 
particles : 
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In  (4),  j (T)  is  the  unidirectional  spectral  intensity  at  the 
observing  point  Q,  T  denotes  kinetic  energy  per  nucleon,  and 


We  first  made  least  squares  fits  to  available  1  AU  spectral 

data  [previous  references]  for  protons  over  the  energy  range 
2  6 

10  <  T  <  10  MeV/nucleon  using  a  logarithmic  version  of  equation 

(4)  in  order  to  find  Kp,  7p,  and  <p^.  Only  <p^  was  allowed  to  be  a 
function  of  epoch.  Taking  pQ  -  <Pp/2  at  each  epoch,  we  then  made 
a  corresponding  treatment  of  the  1  AU  spectral  data  for  alpha 
particles  to  obtain  Ka  and  7  .  The  values  of  the  four  time- 
independent  parameters  in  (4),  as  so  determined,  are: 

Kp  -  1.99  x  105,  7p  -  2.71,  Ka  -  4.27  X  104,  and  7p  -  2.84  when 
j (T)  is  measured  in  (cm  sec  sr  MeV/n)  . 

Then  adopting  these  values  of  Kp,  KQ,  7p,  and  7a,  we  found 
<Pp  and  tpa  separately  for  17  sets  of  proton  and  alpha  particle 
spectra  at  various  epochs  and  radial  distances  as  observed  by 
McDonald  and  collaborators  with  particle  spectrometers  on 
Pioneers  10  and  11,  Helios  1  and  2,  ISEE  3,  and  Voyager  1  and 
reported  in  unpublished  memoranda  [McDonald,  1982]  [McDonald  et 
al . ,  1984],  The  portions  of  the  spectra  below  our  integration 
threshold,  80  MeV/nucleon,  (i.e.,  the  anomalous  component)  were 
ignored . 
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The  most  directly  useful  sets  of  the  spectral  data  of 
McDonald  et  al.  are  the  six  sets  obtained  with  Pioneer  10  because 
these  observations  were  made  on  the  same  spacecraft  as  were  ours. 
Numerical  data  for  these  six  cases  are  given  in  Table  1.  The 
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A.  comprehensive  body  of  omnidirectional  integral 
intensities  of  cosmic  rays  for  T  >  80  MeV/nucleon  is  provided  by 
the  University  of  Iowa  detectors  on  Pioneers  10  and  11  for  the 
period  1972-85  [Van  Allen  and  Randall,  1985].  A  condensed  example 
of  these  data  for  a  matched  pair  of  detectors  (detectors  C)  on 
these  two  spacecraft  is  given  in  Figure  1.  The  two  principal 
plots  in  this  figure  show,  on  a  semi logarithmic  scale  of 
ordinates,  the  counting  rates  of  these  two  detectors  as  a 
function  of  r.  A  straight  line  connecting  two  corresponding 
points  (i.e.,  after  allowance  for  solar  wind  propagation  delay) 
on  the  Pioneer  10  and  Pioneer  11  curves  establishes  the  apparent 
radial  gradient  of  intensity.  Previously  published  results 
[Van  Allen  and  Randall,  1985]  of  this  procedure  are  reproduced  as 
Figure  3.  The  insert  in  the  lower  right  of  Figure  1  shows  the 
time  dependence  of  the  cosmic  ray  intensity  as  measured  by  IMP- 8 
at  1  AU  (courtesy  S.  M.  Krimigis  and  R.  B.  Decker). 

The  counting  rates  C  of  our  detectors  are  proportional  to 
the  sum  of  the  omnidirectional  intensities  of  protons  and  alpha 
particles.  A  consistent  normalization  of  our  counting  rates  to 
the  integrals  of  the  six  available  Pioneer  10  spectra  of  McDonald 


et  al.  (at  various  epochs)  is  obtained  [Table  1}  by  taking  our 

detectors  C  on  Pioneers  10  and  11  to  have  an  omnidirectional 

2 

geometric  factor  g  of  0.157  cm  ,  about  29%  greater  than  our 
previous  best  estimate  of  0.122  cm^  [Van  Allen  et  al.,  1980]. 

This  discrepancy  is  not  understood  but  it  is  not  unreasonable 
because  it  encompasses  the  uncertainties  of  absolute  calibrations 
of  both  Che  GSFC  and  Iowa  detectors.  Adopting  the  normalized 
value  of  g,  we  integrated  formulae  (4)  numerically  from  80 
MeV/nucleon  to  infinity  to  find  "predicted"  values  of  our 
counting  rate  C  as  a  function  of  .  An  abridged  summary  is 
given  in  Table  2. 

In  particular,  the  value  of  C  for  <pp  -  <pa  -  0,  namely  1.83 
count/sec,  is  taken  to  be  C(R) ,  i.e.,  the  interstellar  value. 

This  value  is  represented  by  the  horizontal  line  near  the  top  of 
Figure  1. 

Observed  values  of  C  as  a  function  of  time  and  radial 
distance  r  were  taken  from  our  previous  work  [Van  Allen  and 
Randall,  1985]  and  updates  thereof.  Further  in  accordance  with 
that  work,  we  next  assumed  that  the  radius  R  of  the  modulation 
region  is  given  by 


where  G  is  the  observed  mean  integral  radial  gradient,  having  the 
value  0.0207  (AU)"*^.  The  graphical  significance  of  (5)  is  evident 
in  Figure  1.  We  adopt  C(R)  -  1.83  count/sec  per  Table  2.  Values 
of  R  as  thus  computed  are  plotted  in  Figure  4  for  the  recent 
solar  cycle. 

The  inferred  values  of  R  vary  with  time  systematically  and 
range  from  42  AU  at  the  time  of  minimum  solar  activity  (  ~  1976) 
to  88  AU  about  1.5  yr  following  the  time  of  maximum  solar 
activity  (  ~  1980) .  Also  shown  in  Figure  4  is  a  plot  of  monthly 
mean  sunspot  numbers  [Solar-Geophysical  Data,  1972-1986).  The 
general  correlation  between  R  and  the  sunspot  number  is  striking, 
though  at  the  greater  radial  distances  R  lags  the  sunspot  number 
by  about  1.5  year  [cf.  Van  Allen  and  Randall,  1985). 

Figure  4  presents  the  principal  results  of  this  paper. 


We  consider  that  a  loose  relationship  to  convection- 
diffusion  theory,  as  sketched  in  the  introduction,  is  a  desirable 
aspect  of  our  exposition.  But  such  a  relationship  is  not 
essential  to  our  analysis. 

Indeed,  the  analysis  can  be  framed  in  purely  empirical 
terms  as  follows : 

(a)  Each  of  formulae  (4)  can  be  regarded  as  an  empirical  one 
with  one  (positive)  parameter  <p ,  whose  value  is  a  function  of 
radial  distance  and  phase  of  the  solar  activity  cycle  and  is  zero 
at  the  modulation  boundary. 

(b)  By  fitting  formulae  (4)  to  observed  differential  spectra  of 
protons  and  alpha  particles  at  1  AU  and  at  various  epochs,  with 
the  assumption  that  <p  -  2<pQ,  values  of  Kp,  KQ,  7p,  and  ya  are 
found.  (The  overall  results  are  insensitive  to  the  assumption 
that  <Pp  -  2 <pa  and  are  only  slightly  changed  if  the  alpha  particle 
component  is  ignored.) 

(c)  With  the  time- independent  values  of  Kp,  Ka,  7p,  and  7Q 
from  step  (b) ,  values  of  <pp  and  <pa  are  found  for  a  variety  of 
independently  observed  spectra  at  various  epochs  and  radial  distances 


(d)  As  <pp  and  cpa  tend  toward  zero,  formulae  (4)  reproduce  the 
high  energy  portions  of  observed  differential  spectra  as  repre¬ 
sented  by  formula  (3). 

(e)  Using  values  of  the  spectral  parameters  for  six  Pioneer  10 
cases,  we  find  that  the  ratios  of  the  integrals  of  formulae  (4) 
to  our  observed  integral  omnidirectional  intensities  C  are 
essentially  constant  [Table  1]  over  a  range  of  a  factor  of  two 
in  C. 

(f)  We  then  use  the  corresponding  value  of  C(R)  [Table  2]  to 
infer  values  of  R  by  way  of  equation  (5),  adopting  the  previously 
determined  mean  value  of  the  radial  gradient  G. 

(g)  An  essential  ad  hoc  assumption  is  that  formula  (3)  is  a 
valid  representation  of  the  interstellar  spectra  of  protons  and 
alpha  particles  down  to  T  -  80  MeV/nucleon.  If  this  assumption 
is  wrong  then  the  horizontal  line  in  the  upper  portion  of 
Figure  1  must  be  either  raised  or  lowered  as  appropriate  and  our 
inferred  values  of  R  will  be  increased  or  decreased  in  a 
corresponding  way. 

Our  use  of  a  single  mean  value  for  the  radial  gradient  is, 
of  course,  subject  to  criticism.  It  is  clear  from  Figure  3  that 
piecewise  fits  to  short  portions  of  the  body  of  data  shown  there 
yield  a  substantial  spread  in  values  of  G,  including  negative 
values.  Also,  a  preliminary  analysis  of  Pioneer  10/11  data  for 
1984-1986  yields  G  =  0.013  AU”^.  (As  with  our  previous 


determinations  of  G,  this  analysis  assumes  cylindrical  symmetry 
around  the  polar  axis  of  the  ecliptic  plane.)  The  physical 
reason  or  reasons  for  short-term  fluctuation  in  the  apparent 
radial  gradient  are  not  known  and  we  have  chosen  to  use  the  mean 
value  in  equation  (5)  as  appropriate  to  this  paper,  whose  intent 
is  heuristic  and  not  definitive.  The  consequences  of  using 
different  values  of  G  can  be  visualized  with  the  help  of  time- 
labeled  plots  such  as  those  in  Figure  1 . 

The  smooth  curve  in  the  lower  portion  of  Figure  4  shows  the 
radial  distance  of  Pioneer  10  as  a  function  of  time.  A  specific, 
testable  prediction  of  our  work  is  that  Pioneer  10  will  reach  the 
modulation  boundary  in  1988  and  will  remain  in  its  vicinity  for 
several  years.  A  secondary  prediction  is  that  C  will  increase 
from  its  early  1986  value  of  1.0  to  1.83  count/sec  in  1988. 
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Figure  1.  The  two  principal  curves  give  the  counting  rates  of 
detectors  C  on  Pioneers  10  and  II  as  a  function  of  radial 
distance,  with  the  phase  of  the  solar  activity  cycle  as  the 
implicit  variable.  The  two  slanting  straight  lines  that 
bound  these  curves  are  drawn  for  a  radial  gradient  of 
0.0207  (AU)  ^ .  The  horizontal  line  near  the  top  of  the 
figure  at  1.83  count/sec  represents  our  inferred  value  of 
the  interstellar  counting  rate.  In  the  lower  right  of  the 
figure  is  a  plot  of  the  time -dependence  of  the  cosmic  ray 
intensity  as  measured  by  IMP-8  at  1  AU  [courtesy  of 
S.  M.  Krimigis  and  R.  B.  Decker]. 

Figure  2.  Spectral  data  for  protons  and  alpha  particles  for  the 
second,  fourth,  and  sixth  cases  of  Table  1,  all  from 
Pioneer  10  [McDonald  et  al.].  The  smooth  curves  are  our 
least  squares  fits  of  formulae  (4)  as  described  in  the 
text.  The  two  dashed  curves  represent  asymptotic  spectra 
for  <pp  -  0  and  <pa  -  0,  respectively,  i.e.,  in  interstellar 


space . 
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Figure  3.  A  sample  plot  from  Van  Allen  and  Randall  [1985] 
illustrating  our  determination  of  the  mean  radial 
gradient  of  cosmic  ray  intensity  using  data  from 
Pioneers  10  and  11. 

Figure  4.  The  central  portion  of  the  figure  shows  inferred 
values  of  the  heliocentric  radius  of  the  cosmic  ray  modula 
tion  boundary  in  astronomical  units  (AU)  over  the  recent 
solar  activity  cycle.  The  upper  portion  gives  a  plot  of 
monthly  mean  sunspot  numbers  for  the  same  period.  The 
smooth  curve  in  the  lower  portion  gives  the  heliocentric 
radial  distance  of  Pioneer  10. 
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